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We propose a quantum optical implementation of a class of dissipative spin systems, including the 
XXZ and Ising model, with ultra-cold atoms in optical lattices. Detuned Raman transitions in the 
Lamb-Dicke regime enable engineerable dissipation and a tunable transversal magnetic field. This 
system is an ideal testbed for studying steady-state dynamics of dissipative spin models. For small 
spin chains and weak dissipation, steady-state expectation values show pronounced peaks at certain 
critical system parameters. We show that this effect is related to degeneracies in the Hamiltonian 
and derive a sufficient condition for its occurrence. 



I. INTRODUCTION 

Quantum spin models play a fundamental role for the 
theoretical and experimental study of quantum many- 
body effects. They represent paradigmatic systems ex- 
hibiting, e.g., quantum phase transitions and peculiar 
forms of matter [T] . They also provide toy models for de- 
scription of many solid state systems. Ultra-cold atoms 
in optical lattices [5] have emerged as a system that is 
especially suited to study the low-energy sector of quan- 
tum spin systems with the promise to eventually simulate 
theoretical models in large, controlled quantum systems. 

To observe these effects, coupling to uncontrolled de- 
grees of freedom has to be kept to a minimum, since 
it leads to dissipation and decoherence [31 13] which can 
mask or destroy the quantum effects. But in recent years, 
it has been shown how the coupling to an environment 
can be harnessed to generate useful quantum states [SHI] 
or perform quantum information tasks [SlllOj. Moreover, 
the study of the phase diagram of open quantum systems 
has turned into a fruitful direction itself [TTI - [T6] . 

Our aim in the present work is twofold: In the first part 
of the paper, we propose a scheme to realize a quantum 
spin system using ultra-cold atoms in an optical lattice in 
which both coherent interaction and dissipation can be 
engineered and controlled, enabling the study the non- 
equilibrium and steady-state physics of open and driven 
spin systems. In the second part, we highlight a pecu- 
liar feature of the steady-state diagram for small spin 
chains: in the limit of weak dissipation, abrupt changes 
of steady-state expectation values for certain critical val- 
ues of the system parameters are observed. We explain 
this feature and relate it to degeneracy properties of the 
system Hamiltonian and derive a sufficient condition for 
the occurrence of sharp peaks at critical system parame- 
ters. 



II. PHYSICAL IMPLEMENTATION OF A 
ONE-DIMENSIONAL SPIN CHAIN UNDER 
DISSIPATION 

Ultra-cold bosonic atoms in optical lattices are ideal 
candidates to simulate spin Hamiltonians. Different the- 



oretical and experimental approaches [T7] have been em- 
ployed to simulate quantum spin chains in optical lat- 
tices, for example by optical driving of two hyperfine lev- 
els of cold bosons in the Hubbard regime [15] . Recently, a 
one-dimensional chain of interacting Ising spins has been 
implemented experimentally using a Mott-Insulator of 
spinless bosons in a tilted optical lattice [T^ . 

In the following, we show theoretically how to add engi- 
neered dissipation to the toolbox of these systems [201111] ■ 
Specifically, we show how to implement a system with 
the following properties: [i] dissipative dynamics of Lind- 
blad form, [ii] a tunable magnetic field in a;-direction and 
[iii] an effective spin Hamiltonian such as, e.g., the XXZ, 
Heisenberg or Ising model. In the next Subsections, we 
first introduce the setup and explain qualitatively how 
such a one-dimensional spin chain in a tunable magnetic 
field under engineerable dissipation can be realized with 
cold atoms in optical lattices. In the subsequent Subsec- 
tions we give specific requirements and parameters and 
details of the derivation for [i]-[iii]. 



A. Setup and qualitative description 

The system we consider is an optical lattice populated 
with a single atomic bosonic species. We assume to be 
in the Mott-insulator regime with filling factor 1, where 
the on-site interaction is much larger than the tunnel- 
ing (hopping) between neighboring lattice sites. In this 
regime, the atoms are localized such that each lattice po- 
tential is occupied with one atom. We aim to use the mo- 
tional ground and first excited state of the atom (denoted 
by |0) , |1)[12], respectively) to realize an effective spin- 
1/2 system in each lattice site. To access the motional 
degree of freedom optically, we work in the Lamb-Dicke 
regime where the motion of the atom is restricted to a 
region small compared with the laser wavelength. We 
make use of the anharmonicity of the lattice potential 
and, as explained in the following, of decay of the atoms 
that leads to cooling of the system, to restrict the dy- 
namics to the two-dimensional subspace of {|0) , |1)} [23] 
(see Fig. [T]). For the optical manipulation, we assume 
that the atoms have internal degrees of freedom that can 
be addressed with laser fields. We consider a A-scheme 
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FIG. 1. Relevant level structure and coupling/decay terms of 
single atom trapped at a lattice site. Upper part (Left): In- 
ternal levels of the atom: A system \g), \r), |e), off-resonantly 
driven by lasers. (Right): Motional states in the lattice po- 
tential. Lower part: After adiabatic elimination of |e), an 
effective two-level system with tunable decay rate F and de- 
phasing rate 7 is obtained. 



with two ground states \g) and |r) (both trapped by the 
same optical lattice potential) and an excited state |e). 
The level scheme of the internal states of the atoms is 
shown in Fig. [l] OfF-resonant laser fields drive transi- 
tions between the two ground states \g) and \r) and the 
excited state |e). The system decays fast into the ground 
states, and, as we show below, effectively decays into the 
state \g). Therefore, the atoms are optically pumped to 
the state \g) ® |0) and the states \r) and |e) can be adi- 
abatically eliminated. Eliminating the excited state |e) 
leads to the effective two-level system in the lower part of 
Fig. [1] with designable decay rates. Further elimination 
of the state \r) leads to an effective description in the 
internal ground state \g). The optical couplings by laser 
fields give rise to effective Hamiltonians and effective dis- 
sipation (cooling) in the ground state \g) at each lattice 
site. Details are given in Section |II B[ In summary, we 
obtain an effective two-level system at each lattice site 
with Hilbert space spanned by \g) ® |0) and \g) ® |1) as 
depicted in Fig. [3j 

In the following Sections, we show that engineering the 
optical couplings as above leads to an effective master 
equation for the two-level system |0), |1) that [i] describes 
decay from |1) to |0) and [ii] an effective magnetic field 
in x-direction. In the Mott insulator regime, tunnel cou- 
plings between neighboring lattice wells can be treated 
as a perturbation, which [iii] leads to an effective spin 
Hamiltonian. The resulting master equation 1241 is given 
by 

Pt ^Y.k^^i'^^kPt^k - {^k'^k^Pt}+) -i[H,Pt\ ■ (1) 
Here, — |l)(0|j, is the operator that excites an atom at 



FIG. 2. Effective two-level system ]5)-|r) in the optical lat- 
tice potential with motional states |0) and Choosing res- 
onance conditions as explained in Section [II B| the atoms are 
selectively excited from \g) ® \1) to the state \r) ® |0) and 
spontaneously decay into \g) ® |0). 



lattice site k from the motional state |0) to state |1) and 
o-jT = (o"^)^- The sum runs over all N sites of the optical 
lattice potential. The first part in Eqn. ([I]) describes 
decay from state |1) i nto s tate |0) as depicted in Fig. |3] 
It is derived in Sect ion [TlB The decay parameter A" can 
be tuned by changing the Rabi frequencies of the lasers 
and the dctunings and is given by Eqn. ( |18[ ) in Section 

Hb 



II E The Hamiltonian is given by H 
Hb describes the magnetic field in x-direction given by 



-ffspin, where 



HB=Y.B^{a+ + a-) 



(2) 



where Bx is proportional to an effective magnetic field in 
x-direction. It is derived in Section III CI The Hamilto- 
nian -ffspin describes the spin Hamiltonian 

^^spin = ^ai(CTfeCr^+i + cr^Crf+l) + "20-^(7^+1, (3) 



as derived in Section |II D| The parameters ai and a2 
depend on the properties of the optical lattice potential 
and can be tuned. Therefore, the Hamiltonian -ffspin de- 
scribes the XXZ model, the Ising model or the Ileisen- 
berg model. In the following three Sections, we employ a 
perturbative approach to derive a master equation com- 
prising dissipation of Lindblad form [i] as in Eqn. ([T]), a 
magnetic field in x-direction [ii] as in Eqn. ([2| and an ef- 
fective spin Hamiltonian [iii] as in Eqn. ([s]). For the sake 
of clarity, we derive [i]-[iii] in three separate steps employ- 
ing the approximation of independent rates of variation 
as explained in [25) . 



B. Optical couplings of internal atomic states: 
dissipation of Lindblad form 

In this Section, we show that optically addressing the 
atoms with suitably tuned lasers allows to engineer decay 
as m Eqn. ^. 
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FIG. 3. Decay of the effective two-level system |0), |1) de- 
scribed by the effective master equation derived in Section 
II B The dissipation strength A~ is given in Eqn. (|7|. 



We consider the internal levels \g), |r), |e) of an atom 
at site k. The ground states \g) and \r) can be coupled 
via the excited state |e) by a detuned Raman transition 
of two standing wave laser fields with Rabi frequencies 
fli and fl2. Eliminating the excited state |e) leads to an 
effective coupling between \g) and \r) (see Fig. [T]) with 
rJeff = ^1^2/ Sre where Sre is the detuning with respect 
to |e) (for details see Appendix [A]) . As shown in Fig. [l] 
|r) and |e) are coupled by another laser field (denoted by 
a red arrow). Adiabatic elimination of the excited state 
|e) leads to an effective two-level system (as shown in the 
lower part of Fig. [T]) with states \r) and \g) which has 
designable decay rates F and 7 as derived in ^26] (see 
also Appendix [A| . Thereby, the excited state |e) that is 
broadened by spontaneous emission is eliminated, and 
the effective two- level system allows to resolve the 

motional states |0) and |1) of the lattice potential (note 
that we are in the Lamb-Dicke regime) , as can be seen in 
Fig. [2] Under appropriate resonance conditions that will 
be specified in the following, the atoms are excited from 
state |1) (g) \g) to state |0) ^ \r) and spontaneously decay 
into the state |0) (8) \g) as shown in Fig.[2j Adiabatically 
eliminating the state |r), this corresponds to an effective 
decay from state |1) (g) \g) into |0) ^ \g). Thus the atoms 
effectively remain in the internal ground state \g), such 
that the decay can be written as an effective decay from 
state |1) to |0) as depicted in Fig. [3j 

In Appendix]^ we derive in a perturbative approach 
(that corresponds to an adiabatic elimination of the state 
|r)) a master equation that describes the dynamics of the 
two-level system |0), |1) of the atom. Assuming that the 
driving of level \r) is sufficiently weak such that 

|l]eff|«r,7,^^,l^r|, (4) 

and that the level broadening remains small 

F -f 7 < i., (5) 
the master equation is given by 

k 

+A+ {2a+pta- {a.T^, Pt}+) ^IH^'J, Pt]- (6) 

Here, A"*" determines the strength of the heating terms 
and A~ the strength of the decay terms. For simplic- 
ity, are chosen to be independent of the lattice site k. 



can be made dependent on the lattice site k by choos- 
ing different phases of the driving lasers as explained in 
Appendix |A] Note that A'^ ^ A^ is required for the va- 
lidity of the approximation that restricts to the |0) and 
|1) subspace. A^ and are given by 



(r + 7) 



(F + 7)2 + (J,±i.)2 



(7) 



Here, Sr is the effective detuning given by Eqn. (A5) in 
Appendix |a] r/i = ki/ \j2Mv is the Lamb-Dicke param- 
eter where k\ is the wave number of the laser with Rabi 
frequency fii, M the atomic mass, and v denotes the en- 
ergy difference between the motional state |0) and |1) of 
the lattice potential. The Hamiltonian i/^i'' in the last 



term in Eqn. Q is given by 



(8) 



where Hs describes AC Stark shifts on the motional lev- 
els that are ^ v and are given in more detail in Appendix 
|J| Now, we have everything at hand to implement dissi- 
pation. If 

br w V, 

which can be achieved by choosing the laser frequency 
in br = ujr^ioi accordingly, the strength of the dissipation 
is much larger than the strength of the heating: 



A+ <^A-. 



(9) 



Then, the master equation has only decaying terms and 
is of the form 

Pt =Y.kA-{2<J^Ptat - {a+cjl,pt}+) - I [i/^ff\ptjl.O) 

It describes decay of the atoms from state |1) into |0), 
while the atoms effectively remain in the internal state 
\g). By adiabatic elimination of the internal state |r), we 
have thus shown that a master equation can be derived 
that can be tuned such that it describes almost pure de- 
cay. 



C. Optical couplings of internal atomic states: 
Effective magnetic field in s-direction 

To derive the effective magnetic field in a;-direction, 
we consider a detuned Raman transition. Two standing 
wave laser fields with Rabi frequencies fia and J7f, couple 
the internal ground state \g) and the excited state |e) of 
the atoms as depicted in Fig.|4j The couphng is described 
by the Hamiltonian 



iiah= I]fe^QCOs(fcQXfc) leXslfc 
+ Jib sin (/cfcXfc) |eX.g|j, + h.c. 



(11) 



where fc^, fcf, denote the wave numbers of the lasers and 
Xfc the displacement from the equilibrium position of the 
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FIG. 4. Level scheme and transitions used to implement the 
transverse magnetic field. Left: A detuned Raman transition 
couples the internal ground state \g) and the excited state 
je) of the atom. Right: Adiabatic elimination of the excited 
state |e) leads to an effective magnetic field in a;-direction (see 
Section II C I which drives transitions between the motional 



states |0) and |1) 




FIG. 5. Tunneling between neighboring lattice wells with tun- 
nel amplitudes to and ti. States with two atoms per lattice 
well are treated in perturbation theory in Section [lID[ as the 
on-site interaction is much larger than the tunneling ampli- 
tudes. 



atom at lattice site k. As we are in the Lamb-Dicke 
regime, sm{ki,Xk) ~ Vbif^k ^'^k) ?^ . ^'^^ cos {kaXk) ~ 1. 



D. Effective spin Hamiltonian 

In the Mott-Insulator regime, bosonic atoms trapped 
by a lattice potential with two motional states are de- 
scribed by the two-band Bose-Hubbard model (see 
Appendix [b|) . We denote the on-site interaction by Uqi, 
Uoo and Un |29) and by tg (ti) the amplitudes for atoms 
in state |0) (|1)) to tunnel to neighboring lattice sites. 
We assume that the on-site interaction Uqi, Uoo, Un » 
to, ti such that tunneling between neighboring wells that 
leads to states with two atoms in one lattice well can 
be treated as a perturbation (see Fig. [5]) . Using second 
order perturbation theory 125^ (for a detailed derivation 
see Appendix |b|, we derive an effective spin Hamiltonian 
ffspin given by: 



H. 



(3) 
eft 



-^spin ^" ^ z 



k ' 



(14) 



with 



i^spin = ^ai(a^(T^+i + <jM+i) + a2fT^(T^+i, (15) 



where ai = -2toh/Uoi, a2 = (tg + t?)/(2C/oi) -tg/C/o 
ti/Uii, the magnetic field in z-direction 



ti/Un ~ tl/Uo, 



(16) 



|0Xi| 



and the Pauli spin matrices a^, with ct? 
|lXO|j,. The Hamiltonian given by Eqn. (14| is an effec- 



Under the condition 



five spin Hamiltonian that is tunable by changing the 
lattice properties. If ai, a2 > 0, i/^jj corresponds to 
the XXZ model with a magnetic field in z-direction. If 
the lattice properties can be tuned such that one of the 
tunneling constants to or ti — >■ 0, i?spin is an Ising Hamil- 
tonian with a magnetic field in z-direction. For ai = a2, 
Hspin corresponds to the Heisenberg model. 



\na\,\nb\ « \5, 



(12) 



where 5^ is the detuning of the driving lasers as depicted 
in Fig. |4j the excited state |e) can be adiabatically elim- 
inated and we get an effective Hamiltonian 



(13) 



that describes a tunable magnetic field in x-direction, 
where Ex is proportional to the effective magnetic field 
strength in x-direction is given by 



E. Dissipative one-dimensional spin chain in a 
magnetic field 

In the previous Sections, we showed — for the sake of 
clarity in separate steps — that optical couplings of the 
internal levels can be engineered such that we obtain a 
master equation of Lindblad form [Eqn. ([6])] fulfilling the 
demands [i]-[iii] of tunable dissipation, spin- interaction, 
and tunable transverse field. Combining these results, 
one has to carefully consider the order of magnitude of 
each term. Doing so, we find that the magnetic field in z- 
direction i?^ in Eqn. ( 14) and the Stark shifts in Eqn. (Isl) 



can be of the same order of magnitude as v. Stark shifts 
and Bz lead to an effective energy difference between the 
motional states |0) and |1) given by 



Thus, we have derived an effective magnetic field in 
x-direction that drives transitions between the motional 
states |0) and |1) (as depicted in Fig. [4] on the right), 
while the atoms remain in the internal ground state \g). 



B, 



'+1 



where B^ is defined in Eqn. (15) and s_, s_|_ are the AC 
Stark shifts in Eqn. (|8| (see Appendix |A). Therefore, 
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= 2 ■ 10"^ 
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= 10--^ 
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= 10-" 




Bx/oi2 



FIG. 6. XXZ model with 4 



\a2 and open 



spins, ai — 

boundary conditions under local dissipation of form given by 
Eqn. (201. Upper part: Steady-state expectation value (J^) 



plotted versus the magnetic field B^ja^. Peaks are observed 
that narrow for decreasing the dissipation strength. Lower 
part: Spectrum of the XXZ chain in the magnetic field 
plotted versus B^ja^. Peaks in the steady-state expectation 
value (upper part) appear at crossing points of the Hamilto- 
nian that are marked with black circles. 



with both X and z components. However, only is fully 
tunable, while is large (compared to B^ , ) and re- 
quired to be so by the conditions for adiabatic elimina- 
tion, cf. Eqn. ([4]). However, an effective dissipative XXZ 
chain without any field in z-direction would be advanta- 
geous for observing critical behavior in the steady-state 
dynamics that we study in the next Sections. There- 
fore, we transform to a frame rotating with v. In the 
rotating frame, Bq becomes time-dependent. To obtain 
a time-independent field in cc-direction, the detuned Ra- 
man lasers that lead to the effective magnetic field B^ 
have to be chosen time-dependent, adapted to the rotat- 
ing frame (i.e., suitably detuned from the two-photon res- 
onance) . This then yields a time-independent transversal 
magnetic field, and the master equation in the rotating 
frame is then given by 

-i[ff,pin+i?S,Pt]. (20) 

It corresponds to the master equation given by Eqn. ([T]). 
In summary we have shown how to implement a one- 
dimensional spin chain with nearest-neighbor interaction 
described by the XXZ or the Ising model and a tunable 
effective magnetic field in x-direction under dissipation. 
This system is an ideal testbed for studying steady-state 
dynamics of dissipative spin models as discussed in the 
next Section. Note, that since we are in a rotating frame, 
observables other than the collective spin operator (J^) 
become explicitly time-dependent. 



F. Steady-state behavior: Discontinuous 
steady-state behavior related to spectrum of 
Hamiltonian 



combining all results, the laser detuning (5,. that enters in 



the ^ has to be adjusted to 8^ such that br - 



which means that (5^ = (5r ± (-Bz + s_ — s-|_). 

Then, combining the results from Eqns. ([6]), (131 and 



(14), the master equation reads 

Pt = Efe ^+(2o-A^Pto-fc - {(^l^t. Pt}+) 
+ A-(2a-p(a+ - {a+a-,pj+) (17) 

where the rates are modified by the renormalized (5^: 

(r + 7) 



(18) 



(r + 7)2 + {b,. ± vY 

The Hamiltonian part of the master equation is given by 

i^-i^spi„ + i^B + ;>Efc|iXi|fc, (19) 



where i?spin is given by Eqn. (|T5|) and Hq by Eqn 



The magnetic field in z-direction and Stark shifts have 
been included in v. For (5^ ~ as shown before, decay 
dominates over heating: ^ A^ . Then, the mas- 



ter equation has only decaying terms and Eqn. (17 1 de- 



A particular important characterization of dissipative 
dynamics is through their steady state: if it is unique 
(or distinguished by some conserved quantity) it allows 
for robust preparation of these states. Abrupt changes in 
the steady state as system parameters are varied may sig- 
nal dissipative quantum phase transitions [T2Hl6l [30 l [3T ] . 
We study the steady-state behavior of short spin chains 
under dissipation in a magnetic field in x-direction with 
numerical simulations. We find that the one-dimensional 



scribes a dissipative XXZ spin chain in a magnetic field 



XXZ model with 4 spins as given by Eqn. ( 15 1, where we 
chose as a typical example a\ — \(X2, shows a surpris- 
ing behavior: Changing the external magnetic field in 
x-direction, peaks occur in the steady-state expectation 
values of the collective spin operators J^l^ = J2k '^iJ^ 
for weak dissipation, see Fig. [6] Here, we have considered 
dissipation as in Eqn. ( 20 ) with equal dissipation strength 
on each spin. We find that decreasing the strength of 
the dissipation the peaks become more narrow and each 
peaks height approaches a finite value. For small 7 we 
observe very narrow peaks. This indicates a discontinuity 
in the steady-state expectation values of the spin oper- 
ators. We find, that these narrow peaks appear exactly 
at points where the Hamiltonian becomes degenerate. In 
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FIG. 7. Ising model with 4 spins with open boundary condi- 
tions under local dissipation of form given by Eqn. ( |23[ ). Up- 
per part: Steady-state expectation value (J^) plotted versus 
the magnetic field Bx/o-z. Peaks are observed that narrow for 
decreasing the dissipation strength. Lower part: Spectrum of 
the Ising Hamiltonian in a magnetic field plotted versus 
Bx/ct-i- Peaks in the steady-state expectation value (upper 
part) appear at crossing points of the Hamiltonian that are 
marked with black circles. 



FIG. 8. Ising model with 6 spins with periodic boundary 
conditions under collective dissipation of the form given by 
Eqn. ( 24 1 in the translation and reflection symmetric subspace 
T = R = 1. Upper part: Steady-state expectation value 
(J^) plotted versus the magnetic field Bx/ctz- Lower part: 
Spectrum of the Ising Hamiltonian in the magnetic field Bx 
plotted versus Bx/ctz. 



Numerical studies of discontinuous behavior in 
the steady state 



the following Section we study this phenomenon in more 
generality. 



III. DISCONTINUITIES IN THE 
STEADY-STATE DYNAMICS OF A GENERAL 
CLASS OF ONE-DIMENSIONAL SPIN MODELS 
UNDER DISSIPATION 

In the previous Section we saw that for the one- 
dimensional XXZ model, peaks in the steady-state expec- 
tation values of the collective spin operators appear, that 
are closely related to the spectrum of the Hamiltonian. 
In the following, we study in more generality, indepen- 
dent of a physical implementation, local one-dimensional 
spin Hamiltonians under dissipation of different kinds. 
We present a condition that elucidates the discontinuous 
behavior of the steady state at degeneracy points of the 
Hamiltonian. Then, we study and explain this condition 
in more detail for Ising Hamiltonians. 



We numerically simulate short spin chains. First, we 
study the one-dimensional Ising model with open bound- 
ary conditions, described by the Hamiltonian 

H^H,, + Hb (21) 

with 

i?.. =a3^a^a|+i, (22) 

k 

and Hb as in Eqn. ([2| subject to local or collective decay 
with Lindblad operators cx cr^T or J2k^k^ respectively. 
The master equation describing the full system with local 
dissipation is given by 

Pt = ^7fc (2o-fe PtCr^ - {0-fcCr^,Pt}-H) -i[H,pt]. (23) 
k 

Changing the magnetic field , we find that for weak 
dissipation the steady-state expectation values of the spin 
operators (J^) and (J^) change abruptly at particular 
values of B^., see Fig. [7] Here, we have considered dissi- 
pation as in Eqn. (|23| with equal dissipation strength on 



7 



each spin, 7^ = 7. Decreasing the strength of the dissi- 
pation, i.e., decreasing 7, the peaks become more narrow 
and their height converges to some finite value, while the 
expectation value vanishes elsewhere. For 7 — >■ 0, we ob- 
serve very narrow peaks, which indicate discontinuities 
in the steady-state expectation values of the spin oper- 
ators. We find, that these narrow peaks appear exactly 
at degeneracy points of the spectrum of the Hamilto- 
nian. I.e., to every peak found at some value of = xq 
for 7 — 7^ 0, at least one pair of degenerate eigenvalues 
Ai,A2 of the local spin Hamiltonian H^z can be found, 
i.e., Ai(a;) — A2(x) aX x — xq. Note that the disconti- 
nuities in the steady state at critical system parameters 
are only observed for 7^0. I.e., the (weak) dissipation 
allows us to gain information about the Hamiltonian's 
properties that is not readily accessible in the case of 
7 0. 

This effect can be observed for different kinds of spin 
Hamiltonians such as, for example, the XXZ model (see 
Fig. [6|, both for periodic and open boundary conditions. 
Moreover, changing the type of dissipation, the observed 
behavior does not change qualitatively. E.g., collective 
dissipation, which describes the dynamics of spins all cou- 
pled to the same bath and leads to the master equation 

Pt = l{2J-ptJ+ - {J+J-,pt}+) - i[H, pt], (24) 

where = J2k ^k-> leads to discontinuous behavior 
in the steady-state expectation values as shown in Fig. [8] 
for the Ising model. Choosing an " inhomogeneous" dis- 
sipation which is of the form of the dissipative part in 
Eqn. (23), where now the strengths of the dissipation 



x10 



7fc are different for each spin, peaks can be observed for 
an even larger class of spin Hamiltonians: For 7^, = 7, 
and H = Hh + Hb, where Hh is the Heisenberg spin 
Hamiltonian, we do not observe any peaks. However, if 
we choose different dissipation strengths 7^ for each spin, 
we find peaks at the degeneracy points of the Hamilto- 
nian, as can be seen in Fig. |9] 



B. General condition for discontinuities in the 
steady state 

Since the Liouvillian depends smoothly on the system 
parameters, the observed discontinuities must be related 
to degeneracies in the spectrum of C. As we shall see, 
in the weak dissipation limit they are directly related to 
degeneracy points of the Hamiltonian. 

We consider a system described by the master equation 




FIG. 9. Heisenberg model with 4 spins with open bou ndary 
conditions under local dissipation as given by Eqn. (231 with 
different dissipation strengths 7fe. Upper panel: steady-state 
expectation value (J^) plotted versus the magnetic field B^- 
Lower part: Spectrum of the Heisenberg Hamiltonian in a 
transverse magnetic field B^, plotted versus B^/a^. 



Hq(x) is non-degenerate for x xq. The term Ci con- 
tains dissipative terms and is independent of x. We are 
interested in the limit of weak dissipation 7 — >■ and in 
the change of the steady state at the degeneracy point 

X = Xq. 

The steady state Pssix) is determined by C{x)pss{x) = 
and can be determined perturbatively. The kernel of 
Co{x) is highly degenerate, being spanned by all eigen- 
projectors |Ai(a;)XAi(a;)| of the (non-degenerate) Hq{x). 
This degeneracy is lifted by £1 and the steady state for 
7 — >■ is for X ^ Xq given by 



P^(x)£iP^(x)p3,(x) = 0, 



(26) 



where 



where 



p{t)^Cp=[CQ{x)+^Ci]p{t), 



l^o{x){p) = -~i[H{x),p], 



(25) 



with a Hamiltonian H{x) depending (analytically) on a 
parameter x. For simplicity, we consider the case that 



The possibility of discontinuous behavior of Pss{x) at 
X = Xq arises from the enlargement of the kernel of Cq {x) 
at this point: if Ai and \j become degenerate at x — 
Xq then coherences between degenerate eigenvectors (i.e., 
|Ai(a;))XAj(a;)| ^ j) become stationary at a; = xq. We 
denote by the projector on these additional elements 
in the kernel of £0(2^0) [52]. As we show in Appendix [C| 
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a discontinuity pss{xo) 7^ Imix^xa Pss{x) arises if 

P^£i lim pssix) ^ 0, (28) 

which means that Ci couples the steady state to the 
newly available subspace in the kernel of Cq. For 
simplicity, we made the assumption that the Hamilto- 
nian is non-degenerate for x ^ xq. If the Hamiltonian 
does have degeneracies outside xq, but additional eigen- 
vectors become degenerate at x = xq the argumentation 
follows identical lines, as also in this case, Ci can couple 
the steady state to a newly available subspace P^. 

In Figs. [6] - [9] a large feature appears in the steady- 
state expectation value (J^) around — 0. It narrows 
for decreasing 7, but is not a sharp peak for any of the 
parameters used for 7. This broad peak represents the 
effect of one (or several, cf. Figs. [6] and [9| unresolved de- 
generacies around = 0: Note that for all spin models 
considered, the degeneracy of their respective Hamilto- 
nian is very high at B^ = and is lifted slowly (certain 
eigenvalues touch and do not cross as B^ 0). There- 
fore, the finite values of 7 used in the plots are not much 
smaller than all energy differences and we are not in the 
weak dissipation limit. As 7 is reduced, additional peaks 
are resolved (cf. Fig. |9]). Note also, that some degener- 
acy points do not lead to discernible peaks, e.g. in Fig.|6) 
simply because the corresponding peaks are too small. 



C. Steady-state behavior for Ising Hamiltonians 

To get a better insight into how the condition given 
by Eqn. (28) explains the peaks seen in Fig. [Sj we now 
specialize to the Ising model under collective dissipation 
given by Eqn. (24 1. Then we see that the steady state 



apart from the degeneracy points and the condition for 
discontinuity become very simple. For a detailed deriva- 
tion of what follows, see Appendix [C| 



The Hamiltonian in Eqn. (21 ) with periodic boundary 



conditions is in general degenerate due to translational 
and reflection symmetry. To obtain a non-degenerate H, 
we restrict our consideration to a specific subspace with 
eigenvalue 1 for the translation operator T and the reflec- 
tion operator R [33]. Note that the Hamiltonian is also 
symmetric under the spin- flip operation F = cr^^ , i.e., 
FHF^ = H . Using the properties of £1 and P-invariance 
of iJ, we find that if the system has a unique steady state 
of C{x), it is, in the limit of weak dissipation, given by the 
maximally mixed state cx 1: plugging 1 into Eqn. (26) 
we obtain P^£i(l) = Y.Mx)l\i{x)\ J' \\{x)l\i(^\ 
and flip invariance of H implies (Ai(a;)| \\i{x)) = for 
the eigenstates of a non-degenerate Hamiltonian H (see 
Appendix [C]) . 

Thus if the steady state is unique, it is always maxi- 
mally mixed outside degeneracy points and we see a dis- 
continuity aX X ^ xq if for the degenerate eigenstates 
|Ai(a;o)) , |A2(xo)) we have 



This can be checked to hold for the points at which peaks 
are observed in Fig. [Sj 

For the Ising model in a transverse magnetic field, for 
larger TV the peaks decrease in height, and disappear in 
the limit N — > 00. The spectrum for the Ising model in 
a transverse field is known analytically [33]. For large iV, 
the spectrum is very dense and degeneracy points are so 
closely spaced that peaks are no longer resolvable (and 
vanish in the thermodynamic limit as bands develop). 
Nevertheless, for small spin systems, these features pro- 
vide a method to dissipatively study degeneracies of the 
applicable Hamiltonian - anywhere in the spectrum, not 
just in the ground state. To the extent that £1 is tunable, 
it even provides access to the nature of the degenerate 
states via Eqn. (28). 



IV. CONCLUSIONS 

We have shown that using cold atoms in an optical lat- 
tice in the Mott-insulator regime, dissipative spin chains 
with Hamiltonians such as the XXZ model, the Ising 
model or the Heisenberg model can be realized. Optical 
driving of internal atomic states allows for the realization 
of a tunable transversal magnetic field and engineered 
dissipation. 

This system is an ideal testbed for studying steady- 
state dynamics of dissipative spin models. We have dis- 
covered a peculiar feature of the steady state diagram for 
small spin chains: in the limit of weak dissipation, the 
expectation values of the collective spin operators exhibit 
abrupt changes that hint at discontinuities in the steady 
state. These discontinuities occur at degeneracy points 
of the Hamiltonian. We have studied this phenomenon 
for different spin models with open and periodic bound- 
ary conditions subject to individual and collective dissi- 
pation. Finally, we have presented conditions that elu- 
cidate the discontinuous behavior of the steady state at 
degeneracy points of the Hamiltonian. Therefore, mea- 
surements of the steady state dynamics of cold atoms in 
optical lattices would allow to draw conclusions on the 
spectrum of the respective spin model. 



Appendix A: Derivation of effective dissipative 
master equation 

The internal levels of the atom that we consider are 
Ig), |r), |e). Adiabatically eliminating the excited state 
|e) as discussed in Section II B we get an effective two- 
level system \g) and \r) that is coupled with the effective 



(All J^|A2> 7^0. 



(29) 



Rabi frequency fJcS as depicted in Fig. [TJ 

In the following, we derive in detail the master equa- 
tion given by Eqn. ^ in Section II B[ The internal levels 
of the atom that we consider are \g) , |r), |e), as depicted 
in the upper part of Fig.[T] The states are coupled 

by a detuned Raman transition via the excited state |e) 
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by two standing wave laser fields. The coupling is de- 
scribed by the Hamiltonians 



Hii = ^ f^i cos (kiXk){\e){g\f. + h.c), 

k 



(Al) 



and 



Hi2 = 51 ^2 Sin (fc2Xfc)(|rX5|fc + h.c), (A2) 



where ili and VL2 are the Rabi frequencies of the two 
lasers and ki and k2 are the wave numbers of the lasers 
and k denotes the lattice site. Xk is the displacement 
from the equilibrium position a;^ of the atom at lat- 
tice site k. The phase of the lasers is for simplic- 
ity chosen such that cos[ki{xk + x'^)] = cos{kiXk) and 
cos [k2{xk + a;°)] = sin(fc2a;fc). Choosing different phases 
of the lasers makes in Eqn. ([7| dependent on the lat- 
tice site k. Adiabatic elimination of the excited state |e) 
leads to an effective coupling 



Hi=Y, ricftr/iK + 4)(kX5l + h.c), 



(A3) 



with f^ofi = ^1^12/ 5re where J^e is the detuning with re- 
spect to |e) and rji is the Lamb-Dicke parameter. Here, 
we have expressed the deviation from equilibrium po- 
sition, Xk, in terms of harmonic oscillator operators 
truncated to the two lowest lying levels sm{kiXk) ~ 
Vii^^k + <^k) where cr+ = |lXO|j^ and a- = \OXl\^ and 
cos(fc2a;fc) « 1. The effective coupling with Rabi fre- 
quency r^cfi between states |r) and \g) is shown in Fig. [l] 
Coupling the state |r) to the excited state |e) with a 
third standing wave laser field with Rabi frequency flen 
depicted with a red arrow in Fig. [T] we can derive an 
effective two-level system with designable decay 

rates as done in Here, we briefly review this result. 
Following the upper level |e) can be adiabatically 
eliminated if the saturation parameter for the transition 
|r) and |e) is small 



52 

'-'re 



(r,, + reg)2/4 



< 1. 



(A4) 



According to [35] , the effective detuning and the effective 
decay rates are given by: 



Sr = S, 



gr 



^[(r,, + re.)/2]2 + 52^' 



^ eg H" ^ er 



[(r,, + re.)/2]2 + 52^ 



(A5) 



(A6) 



(A7) 



see also the lower part of Fig. [T] The effective two-level 
system |(7)-|r) with the effective decay rates F, 7 and the 



effective detuning 5r is the starting point of the following 
discussion. The full Hamiltonian describing the system 
is given by 



Hi+Ho, 



(A8) 



where Hi describes the atom-light interaction given by 
Eqn. (A3 1 and Hq defines the energies of the system 



Ho=j:k^r\r){r\, + ,y\m,. 



(A9) 



The effective dynamics of the system can be derived 
considering contributions to the Liouvillian up to second 
order in a perturbative approach. The full system is de- 
scribed by a Liouvillian given by: 

p{t) = i^o + Ci)p{t), (AlO) 

where £0 is given by 

Cop{t) = r (2 IgXrIfe Pit) \rXg\, ~ {|r)(r|, , p{t)} 
+ l(2\r){r\,pit)\rXr\,-{\r){r\,,p{t)}^ 



i[Ho,p{t)]. 



(All) 



The first part of the Liouvillian is the decay part with 
the effective decay rate F from state \r) to \g) and the 
dephasing rate 7. The projector 



^^.^IffXffl^doXoi + liXil) 



(A12) 



is stationary under Cq. The perturbative part of the 
Liouvillian is given by 



Cip{t) = -t [Hi,p{t)] 



(A13) 



where Hi is given by Eqn. ( A3 1 and describes the in- 
teraction of the two-level system with the effective laser 
field. Treating Ci as a perturbation, we derive an effec- 
tive Liouvillian in the stationary subspace of Cq. The 
projection onto this subspace reads 



¥p{t) = PCPp{t) + ¥£Qp{t), 



(A14) 



where Fp = \g){g\ ® (|0X0| + |lXf |)p \g){g\ ® (|0X0| + |lXl|) 
and Q = 1 — P. Projecting onto the subspace we want to 
eliminate we get 



QCpit). 



(A15) 



In the following, we integrate Eqn. ( |A15 1 to get the time 
evolution of the density matrix in the fast space. 



We insert the result in Eqn. (A14| to get an equation of 



motion for the density matrix in the slow space. There- 
fore, we first go into the interaction picture, where the 
density matrix is given by p{t) = e-^^^pit). The equa- 
tion of motion in the fast space reads 

Qm = QWi{t)p{t), (AI6) 
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with Wi{t) = e'-'°*£ie'-'°* . Solving this equation by iter- 
ation [35| we get 



= 

■t 



Cot 



dsWi{s)Vp{0) 



dsi / ds2Wiisi)Wi{s2)¥m 

"'0 



(A17) 



At time t = 0, p{0)—p{0) and we assume that at t = 0, 
the population is in the ground state, i.e., p{0) = Pp(0). 
Higher order integrals are neglected making the assump- 
tion that 



Ifieffl < r,7, \Sr\,iy. 



(A18) 



We denote the first integral in Eqn. (A17) by Ri{t) and 
the second integral by i?2(i) such that 

Qpit) = Ri{t)+R2{t). (A19) 

Inserting in Eqn. ( |A14 ) leads to 

Pp{t) = FCFp{t) + ¥£aRi{t) + FCiRi{t) 

+ VCoR2{t)+VCiR2{t). (A20) 

The term PCoRi{t) ^ 0, and PCiR2{t) is a third or- 
der term and can be neglected. Neglecting terms rotat- 
ing with exp(±ii/i) we get the master equation given by 
Eqn. ^ with AC Stark shifts given by 



where 



s± = n^grji 



{5r ± v) 



(r + 7)2 + {5r ± 



Appendix B: Derivation of the spin Hamiltonian 

In the Mott-Insulator regime, bosonic atoms trapped 
by a lattice potential with two motional states are de- 
scribed by the two-band Bose-Hubbard model 



(Bl) 



Here, the sum runs over the N sites k of the optical 
lattice. The unperturbed Hamiltonian is given by 

where Uxx' is the on-site repulsion of two atoms on lat- 
tice site k, where one atom is in motional state \x) and 
the other one in \x') with x,x' = 0, 1, respectively. The 
operator fikx — I^^X^^Ij, counts the number of atoms at 
lattice site k in the motional states x = 0,1 and v is the 
energy difference between ground and first excited mo- 
tional states. We assume the system to be prepared in 
the ground state |0). Due to the anharmonicity of the 



potential, we do not leave the subspace of n = and 
n = I excitations. 

The perturbative part of the Hamiltonian describes the 
tunneling between neighboring lattice sites and is given 
by 



(B2) 



Here, the operators Ckx with a; = 0, 1 are bosonic destruc- 
tion operators for atoms in the two motional states |0) 
and |1) at lattice site k. to(^i) £^re the tunneling ampli- 
tudes from state |0) (|1)) at lattice site k to state |0) (|1)) 
at fc -I- 1. 

As the on-site interaction Uxx' S> to,ti, tunneling be- 
tween neighboring wells that leads to states with two 
atoms in one lattice well can be treated as a perturba- 
tion. For that, we consider two neighboring lattice sites 
k and fc -I- 1 and write the effective Hamiltonian in the 
basis of eigenvectors of Hq, jx^, j/fe+i), where for example 
|0/£, Ifc+i) is the notation for the state with one particle 
in well k in state |0), and one particle in well fc -I- 1 in 
state |1). In perturbation theory [251, the second-order 
effective Hamiltonian can be evaluated in the following 



way: 



{xk,yk+i\H^s \x'k,y'k+i) = 
^ {xk, yk+i\Ht \x) -j^ {x\ Ht \x'k,yk+i) 



2 ^ 

X 



(A21) where 



1 

E' 



E. 



xy 



E. 



x'y' 



E. 



and Ix) are eigenstates of Hq with two particles in 
one well (and no particle in the other one). E^y and 
E^ are the unperturbed energie s {xk , yk+i \ Hq \xk, yk+i), 
(xl Ix) etc. Evaluating Eqn. (B3) leads to the effective 



(3) 

spin Hamiltonian H^^ given by: 



-"off ~ ^spm 



k ' 



(B3) 



with 



fc 



«i(CTfef^fe+i + ^^fef^fc+i) + a2fT|a^+i. (B4) 



Here, 



0.2 



2Uoi Uo t/i ' 
and Bz, the magnetic field in z-direction is 



Un 



Thus, we have derived an effective XXZ-spin Hamiltonian 
with a magnetic field in z-direction. 
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Appendix C: Condition for discontinuous behavior 

Here, we first derive a general condition for tlie dis- 
continuous beliavior in tlie steady state at a degeneracy 
point of a large class of spin Hamiltonians. Then, we fo- 
cus on more specific Hamiltonians. We study the steady 
state of flip-invariant Hamiltonians outside the degener- 
acy point and, starting with the general condition for 
finding discontinuities in the steady state, we derive a 
more precise condition for flip-invariant Hamiltonians. 



1. General condition for discontinuities in steady 
state 

Here, we derive a general condition for discontinuous 
behavior in the steady state at the degeneracy point x — 
xq of a general Hamiltonian H , where H = H{x) is an 
analytic function of x. We consider a system described 
by the master equation 



(Cl) 



where the Hamiltonian part of the Liouvillian is given by 
Co{x) = Co = —i[H{x), ■] and depends on a parameter 
X, and the local decay Liouvillian is 



k 



[2a-p(t)a+-{fT+a-,p(t)}J. (C2) 



First, we want to describe the system outside the degen- 
eracy point, i.e., for x xq. We assume that in the vicin- 
ity of xq, the Hamiltonian is nondegenerate (for x xq) 
and that the dissipation is weak. The steady state Pss{x) 
defined by {Co{x) + Ci)pssix) = is, in the limit 7 — > 
given by 



p^(x)/:iP^(x-)p., = 0, 



(C3) 



where P^{x) is the projector onto the kernel(£o)- As 
the kernel of £0 is spanned by the eigenprojectors 
\Xi{x)){Xi{x)\ of H we have for arbitrary A: 



(C4) 



where \Xi(x)) are eigenstates of the Hamiltonian H{x) 
which is assumed to be nondegenerate. 

Now, let us consider the case that at a: = xq, the 
Hamiltonian has a degeneracy point at which two or more 
eigenvalues cross. At this degeneracy point, we expect an 
discontinuous behavior of the steady state that leads to 
the peaks we observe in our numerical simulation (see 
Figs. [6]|9]). At X = Xq the projector onto the kernel of £0 
has to be extended. It now also projects onto coherences 
between eigenstates of H: |Ai),|A2) which are eigenvec- 
tors to the degenerate eigenvalues Ai = A2. Therefore 
the projector on the coherences reads: 



It is convenient to define a continuous extension of the 
projector P''^ a,t x = xq which reads 



P^{xo) = lim P^{x). 

X^Xq 



(C6) 



Thus, at X = Xq, the full projector onto the kernel of £0 
reads P^(a;o) -I- P^- Now the condition for the steady 
state Pss{x — xo) at the degeneracy point is given by 

[P^(:eo) + P^] £1 [P''(^o) + P^] Pss{xo) = 0. (C7) 

We want to find a sufficient condition for the steady state 
to change discontinuously. This means that 



Pss{xq) - lim pss{x) ^ 0, 

X—>-Xq 



(C8) 



where lim pss{x) is the continuous extension of 

Pss[x)yx ^ a;o to X = Xq. 

A discontinuity in the steady state as described by 
Eqn. ( |C8[ ) can occur only if 

[P^(xo) -t- P^] £1 [P^(xo) + P^] lim p,,{x) ^ 0, 

(C9) 

holds, since otherwise the continuous extension 
lim Pss{x) would be a steady state as well. The 



last part of Eqn. ( C9 1 can be simplified using 



[P^(xo)+P^] lim pss{x)^ lim p,,(a;), 

X^Xo X^Xq 

which holds since lim Pss{x) is per definition in the 

X^Xq 

space onto which P^(a;o) projects and is orthogonal 
to that space. By Eqn. (C3| we then see that Eqn. (C9) 
reduces to the condition 



P^£i lim pss{x) ^ 0. 



(CIO) 



If this condition is fulfilled, then Pss{xq) — 
lim Pss{x) 7^ which means that the steady state shows 

x—¥xo 

discontinuous behavior at the degeneracy point x — xq. 



2. Condition for discontinuous behavior for Ising 
Hamiltonians 



Here, we want to get a better insight how the condition 
given by Eqn. (CIO) relates to the peaks observed in our 



numerical simulation. In the following, we will apply it 
to the Ising model in a transverse magnetic field. In the 
numerical simulation (see Fig. |8| for the Ising Hamilto- 
nian with periodic boundary conditions under collective 
dissipation described by Eqn. (24), we restrict our con- 



P^^ = |AiXAi|A|A2)(A2H-h.c. 



(C5) 



sideration to a specific subspace with eigenvalue 1 for 
the translation operator T and the reflection operator R: 
T = R = 1. First, we want to prove that if the steady 
state is unique, it is the fully mixed state outside the de- 
generacy points as indicated by our numerical simulation. 
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Then, we show that starting from the condition given by 
Eqn. (CIO), specialization to the Ising model allows to 



derive a more precise condition for finding a discontinu- 
ity in the steady state at the degeneracy points. 

First, we show that 1 satisfies {Cq + Ci)l = outside 
the degeneracy point x ^ xq. Therefore, for systems with 
unique steady state, it is given by the fully mixed state 
for a; ^ xo in the limit 7 — >■ 0. The Hamiltonian given 
by Eqn. (21 1 is assumed to be non-degenerate for x xq 
and invariant under the spin flip operator F = crf^, i.e., 
FHF^ = H . Thus, we want to show that for x ^ xq: 



V^Cil = 0, 



(Cll) 



where is given by Eqn. (C4|. Then, 

P^£i(l) = 27(J-J+- J+J-)cx7J^ (C12) 



Therefore, Eqn. (Cll) reads 



P^£i(l) - P^J^ - I^'XA.I r = 0(C13) 



If we can show that Eqn. (C13), 



(A,| J^|A,) =OVi, 



then we have shown that the fully mixed state is a steady 
state of our system outside the degeneracy points of the 
Hamiltonian. As the Hamiltonian is nondegenerate and 
invariant under the fiip operator F , the eigenvectors of H 
are eigenvectors of F: F = at \Xi). Let \a) denote an 
arbitrary eigenvector of H with _F-eigenvalue a. Since the 
spectrum of F is{±l}, we have \a) = \a) = aF\a). 
Moreover, the flip F changes the sign of J^, i.e, and 
F anticommute: {F, J^}+ = 0. Therefore, we can write 



Eqn. (C14) as 



(a| \a) = a(a| J^i^|a) 
= -a{a\FJ^\a)^-{a\J''\a), (C15) 

where we have used that — \. It follows that 



= 0. 



(C16) 



Consequently, P^£i(l) — and we have shown that in 
the limit of weak dissipation, the steady state, if it is 
unique, is the fully mixed state. For the Ising model 
with up to 8 atoms and collective dissipation, we know 
from our numerics that the steady state is unique. 

To see that the steady state shows discontinuous be- 
havior at the degeneracy point x — xq, we need now only 
to show that the fully mixed state is not the steady state 
of the system. Thus, we need to show that 



P^£i(l) 



E 



A,XA,| J^|A,XA,| ^0 (C17) 



where P'^ is given by Eqn. (C5) and sums over the 
labels of degenerate e igenv alues. Since the eigenvectors 
are orthogonal, Eqn. (C17) holds if 3i ^ j such that 



^ ^ „ (C18) 

(C14) Therefore, Eqn. (CIS I gives a condition for finding dis- 



(A.|J^|A,>^0. 



continuous behavior of the steady state of the Ising model 
in a transverse field under collective dissipation. Note 
that this derivation can be easily extended to all non- 
degenerate Hamiltonians that are fiip-invariant. 
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